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Generalization of Gaussian trial wave functions in quantum molecular dynamics 
models is introduced, which allows for long-range correlations characteristic for 
composite nuclear fragments. We demonstrate a significant improvement in the 
description of light fragments with the correlations. Utilizing either type of Gaus- 
sian wave functions, with or without correlations, however, we find that we cannot 
describe fragment formation in a dynamic situation. Composite fragments are 
only produced in simulations if these fragments are present as clusters in the sub- 
structure of original nuclei. The difficulty is traced to the derealization of wave 
functions during emission. Composite fragments are produced abundantly in the 
Gaussian molecular dynamics in the limit ?i — > 0. 



1 Correlated Molecular Dynamics 

This talk is devoted to the description of fragments and their production in 
models of molecular dynamics with and without correlations. The considered 
correlations are long-range, such as characteristic for composite fragments, 
different from the short-range correlations discussed by H. Feldmeier. In the 
talk I shall, in sequence, discuss the introduction of correlations into molecular 
dynamics, describe the corresponding equations of motion, and discuss results 
for light clusters: deuteron and alpha. Subsequently, I shall turn to reactions 
and discuss results from simulating the decay of an excited system formed 
in reactions and results from collision calculations. I will end the talk with 
conclusions. 

In quantum molecular dynamics models, the exact iV-particle wave func- 
tion of a system is replaced by a trial wave function used in the time-dependent 
variational principle. The different molecular approaches include Quantum 
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Molecular Dynamics Antisymmetrized Molecular Dynamics (AMD)l 

Extended Quantum-Molecular Dynamics (EQMD) □, and Fermion Molecular 
Dynamics (FMD) Em. A common feature of these approaches is that the trial 
wave function is a product, either antisymmetrized or not, of the single-particle 
Gaussian wave functions 

ip„(x) — n exp {—a v (x - f v ) 2 + ip v x). (1) 

In effect of the product assumption, correlations between nucleons, other than 
those possibly induced by the antisymmctrization, vanish. The molecular ap- 
proaches, otherwise, differ in some important details. Thus, QMD is the most 
classical of the approaches; the trial wave function is primarily used in the 
derivation. The width in the wave function (|l]) just acts to increase the range 
of effective interaction; the dynamics of the width is suppressed. No antisym- 
metrization is carried out and its effects are simulated using a Pauli potential. 
Within AMD the antisymmetrization is carried out. In EQMD the spreading 
of wave functions with time is taken into account. In FMD both the antisym- 
metrization is carried out and the width dynamics is included. 

In any of the existing approaches, the center-of-mass wave function for 
any nucleon cluster is localized since, for a product wave function, the center 
of mass gets necessarily localized once the internal coordinates are localized. 
When a low-energy nucleon is emitted from a nucleus, it becomes delocalized 
getting rid of its localization .energy. That has been, in fact, observed in the 
FMD simulations of reactionsQ. The derealization, however, cannot take place 
in the molecular dynamics models for clusters, since their intrinsic state must 
stay localized. For light clusters such as deuterons or alphas, the cm localiza- 
tion energy is of the order of AE ~ (10 — 20) MeV. Given that temperatures 
are low when the fragments are expected to be produced in the low-energy 
reactions, T < 10 MeV, the cm localization energy can result in significant 
suppression factors for the fragment emission, e~ AB / T < (1/10 — 1/3). 

Aiming at the description of fragment production in reactions, we consid- 
ered ways of solving the problem of the artificial cm localization caused by 
the assumptions within the molecular dynamics. Let us first take a look at 
a deuteron Gaussian wave function, in which the intrinsic and cm motions are 
separated, 
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The wave function that is a product of Gaussians in the cm and relative co- 
ordinates represents an exponential of the bilinear form in the coordinates 
of individual nucleons in (||). The case of a delocalized deuteron with to- 
tal momentum tending to zero corresponds to A' — > oo in (^) and a small 
finite A. The two nucleons are then correlated: the cm position may be any- 
where in space but if one nucleon is found somewhere, then the other nucleon 
is within the distance of the order of A away from the first. On the other hand, 
the wave function reduces to a product of single-particle wave functions only 
when A = A'. 

Equation (|J) suggests to use a generalized Gaussian as a trial wave function 
for N particles, 

(xi, . . -,x N \^} = AT exp [-Aij (Si - fi) • (xj - Tj) + ip t ■ x,] , (3) 

where A is the width matrix. The generalization (J3[) of the product of Gaussian 
single-particle wave functions allows for a full separation of the cm and intrin- 
sic coordinates for any nucleon cluster and introduces correlations between 
different nucleons, i ^ j. 

The width matrix A in (|J) is a symmetric, complex, N x N matrix with 
a positive-definite real part. The parameters {r*i} and {pi}, respectively, repre- 
sent centroids and momenta for wave packets and the set {A^, fi,fi} is related 
to various expectation values, including those that quantify the correlations, 
by 

(xi) = n (4) 
HV 4 )=k (5) 

<(*,). = (fi) a (r)) b + S -f (ReA)-- 1 (6) 

(i(Vi)a i(Vj)b) = (pi) a (Pj)b + 8 ab {A* {Re Ay 1 A)a , (7) 



where a and b are indices for cartesian coordinates. For standard molecular 
dynamics, the width matrix is diagonal in particle indices 

Aij = A, 5ij . (8) 

2 Equations of Motion 

Equations of motion for the parameters {Aij, f{,Pi} = {qp.} follow from the 
variational principle 

S [ dt(m^--H\V) =0. (9) 



3 



(For unrestricted variation, the principle yields the exact Schrodinger equa- 
tion.) The equations take a form 



Avpip = —=—(H) 



(10) 



where the skew-symmetric matrix A represents the product of wave function 
derivatives 



/ d 

A vli =2 Im( — * 
\ dq v 



_d_ 



(11) 



To solve the equations of motion, the matrix A needs to be inverted. When 
no antisymmetrization is carried out, the structure of the matrix is 
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(12) 



where the rows and columns refer, in sequence, to the parameter subspaces 
{RcA}, {ImA}, {r}, and {p}, I is the identity in iV-particle space, and B is 
a quadratic matrix of the dimension N(N+l)/2. The matrix B can be inverted 
analytically to yield a result in terms of the matrix elements of A. Because of 
this analytic inversion, the final effort in solving the equations of motion scales 
with particle number like ~ TV 4 ; with the antisymmetrization the effort scales 
like ~ N 6 . 

Explicitly, the equations of motion take the form 



— Re A = B 
dt 



d(H) ^ 
dpi 

-i d(H) ^ 
d Im A ' 



Pi = 



9(H) 



d T , _! d(H) 
— 1mA = -B 1 — - — '- 
dt dReA 



(13) 



(14) 



Contribution to the rhs of the equations for A from the kinetic energy can be 
written as 



B 



-i d(T) 



Im A 
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Re A 



dim A m ' dReA 
From ( |l4] ) and (15), the A- matrix equation for free particles is 



A 2i A A 

A = A A 

m 



(15) 



(16) 
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with a solution 

A(t)=A (l + -(t -t o )A ) . (17) 

For initial conditions in a Gaussian form, the solution ( p"7j ) represents an exact 
solution to the Schrodinger equation, with the spreading in time in the rela- 
tive and cm coordinates for any set of nucleons. If the matrix A is initially 
diagonal in these coordinates within the TV-particle space, it will stay diagonal 
in these coordinates later on. When an interaction is included on the rhs of 
equations (|l4|), it affects the evolution of the matrix elements and spreading 
for the intrinsic motion, but not for the cm motion. When, though, the ma- 
trix A is always restricted to a diagonal form in the single-particle coordinates, 
then the intrinsic and cm motions get coupled. If spreading in the intrinsic 
coordinates is limited, so becomes the spreading in the cm coordinates. The 
intrinsic and cm energies are not separately conserved. Consequences of the 
coupling will be illustrated for light clusters. 



3 Light Cluster Description 

Figure [l] shows the evolution of different energies for an isolated deuteron, 
within the standard uncorrelated wave-packet dynamics (u) and within the 
correlated dynamics (c) with nonvanishing off-diagonal matrix elements in A, 
following Eqs. (|l^) and (|l4|). In the correlated dynamics the various energies 
for the deuteron stay constant. In the uncorrelated dynamics, the intrinsic 
motion taps on the energy from the cm motion. The intrinsic energy increases 
with time and at t ~ 12 fm/c the deuteron becomes unbound. Figure |^ 
shows the evolution of the elements of the width matrix for the deuteron in 
the uncorrelated (dashed line) and correlated (remaining lines) dynamics. In 
the correlated dynamics the intrinsic width (corresponding to the dot-dashed 
line) stays constant, while the cm width (solid line) behaves like one for a free 
particle with twice the nucleon mass. The single-particle element An (dotted 
line) drops with time and then stabilizes. In uncorrelated dynamics the single- 
particle element drops continuously; the deuteron dissolves. 

The intrinsic state of an alpha particle initialized in the lowest state of 
intrinsic energy is stationary in the correlated wave-packet dynamicsQ. The cm 
state evolves in the same manner as the state for a particle with four times 
the nucleon mass. In the uncorrelated dynamics, on the other hand, the alpha 
particle, initialized as in the correlated dynamics, pulsates with time due to the 
energy exchange between the intrinsic and cm degrees of freedom. The energy 
acquired from the cm is not large enough in this case to dissolve the particle. 
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Figure 1: The dashed, solid, dash-dotted, and dotted lines show, respectively, the kinetic, 
total, internal, and potential energies, as a function of time, for a free deuteron in corre- 
lated (c) and and uncorrelated (u) dynamics. 



The results for the isolated particles show the benefits of using the corre- 
lated over the uncorrelated wave functions in the fragment description. We now 
turn to fragment production. 



4 Reactions 

Two types of simulations relevant for heavy-ion reactions will be discussed. 
One will be the decay of an excited system formed in the reactions. Second 
type will be complete reaction simulations followed from an initial state. 

In our decay simulations we initialized the system (e.g. of mass A = 80) in 
the state of a local thermal equilibrium with radial flow: The centroids fi were 



G 



randomly selected within a sphere of radius R — (Airp/^)" 1 (R = 5 fm). The 
momenta pi were selected randomly from a hnitc-T Fermi-distribution (T = 
5 — 12 MeV). To every pi, a radial-flow momentum was added, proportional to 
position pi oc r\ (E co i =0 — 25 MeV/nucleon). The width matrix was either 
taken proportional to unity (Ax = 1.9 fm), or its diagonal and off-diagonal 
elements were generated randomly according to a distribution, or, finally, the 
matrix elements were equilibrated by keeping the system, at first, in a narrow 
external oscillator potential. 

The initial (panel (a)) and late stages (panels (b)-(c)) of a system prepared 
in the described manner are shown in Fig. ^|. In the calculations with a dynamic 
width matrix, the excited system emitted a number of nucleons but no clusters, 
neither of intermediate mass (IMF), nor lighter such as deuteron or a particle. 
There are some centroids in Fig. ||(c) that are close together, but the associated 
widths arc so large that the interaction cannot hold the nucleons together and 
the centroids will separate with time. 

When you get a final state of simulation without clusters you think first 
that it is a coincidence and you repeat the simulation with somewhat changed 
initial conditions. Again you find no clusters. Then you do that again and we 
have done that tens of times, maybe close to a hundred, and never saw any 
fragments, even an alpha. We observed, in fact little difference in the evolution 
for centroids between the uncorrelated dynamics and the correlated dynamics 
with the A-matrix initialized in a diagonal form. The simulations initialized 
with finite off-diagonal elements in A also resulted in one residue and many 
emitted nucleons, but no IMFs or alphas. 

As the primary reason for the lack of fragment production in the simula- 
tions we find the spreading of the wave packets with time. By the time the 
system reaches low densities and the fragments are expected to form, the pack- 
ets get so spread out that the mutual interaction between the packets is unable 
to shrink them into clusters. Clearly, in the past, the production of IMFs and 
light clusters has been observed in the QMD calculations. Indeed, if we sup- 
press the width dynamics (corresponding to Ti — » in our equations of motion), 
we begin to observe fragment production. This is illustrated in Fig. ||(b) that 
shows a late stage of the excited system evolved following the dynamics with 
a frozen wave-packet width. Two IMFs as well as two deuterons can be seen in 
the panel. Additional small clusters have been emitted before the shown time 
and left the displayed spatial region. 

Our results on cluster production for the dynamic width matrix may seem 
in contradiction to the FMD results a with even multifragmentation events 
reported in nuclear collisions. Indeed, in our own simulations of collisions 13, 
we observed the production of clusters, but only if these were present in the 
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substructure of original nuclei and did not manage to dissolve in the reaction. 
We never observed clusters that were formed in a reaction. For example, in 
the ground state of 12 C, for our interaction, the centroids form three a-type 
clusters of four nuclcons each. If we collide two 12 C nuclei, one or two of them 
may break up into three a particles. On the other hand, if we collide two 40 Ca 
nuclei, with no internal substructure, one or two residues are formed and many 
nucleons are emitted, but no clusters of intermediate mass. 

5 Conclusions 

In general, for uncorrelated wave functions, the internal state cannot be local- 
ized without localizing the center of mass. Thermal estimates indicate that the 
localization can act to suppress fragment production. In the dynamics, the use 
of uncorrelated wave functions leads to unphysical exchange of energy between 
the cm and intrinsic degrees of freedom. These deficiencies are eliminated in 
the dynamics utilizing correlated wave functions. The correlations that we 
have introduced are of such nature as expected for fragments and they allow 
for the separation of any number of nuclcons in a wave function. The benefits 
of the correlated wave functions have been illustrated with the examples of 
deuteron and a particle. 

Despite the benefits, in reaction simulations with dynamic widths, with or 
without correlations, we observe the cluster production only when the clusters 
are present in the initial state. The absence of new clusters is associated with 
large spreading of the wave function at the time when the reacting system 
expands and the new clusters should form. While there is nothing unphysical 
in the spreading of the wave function, as such, in the simulations, in reality 
the interaction would be capable of pinching portions of two or more packets 
and forming fragments out of the portions of these packets. The Gaussian 
parametrization of the wave function does not permit for the process. 

On the basis of the reaction simulations, we conclude that, in a successful 
description of fragment production, the wave function must have a flexibility 
to change over distances comparable to the interaction range, at a time when 
the fragments are formed. The widths of the packets in AMD or QMD may 
already act to suppress the fragment production. 
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Figure 2: Evolution of the elements of the width matrix for the internal ground state of a free 
deuteron. The dashed lines shows the evolution of the element in the case of an uncorrelated 
wave function. The dash-dotted, solid, and dotted lines show, respectively, the evolution of 
the internal, center-of-mass, and diagonal elements in the case of a correlated wave function. 
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Figure 3: Particle centroids (dots) in the configuration space for A = 80, in the initial 
state (a) and at t = 200 fm/c for the static wave-packet width (b), and for the dynamic 
width matrix (c) and (d). In (c), the initial flow energy is lowered, compared to (b) and (d), 
by the initial energy content in the localization of the wave packets. The circles indicate rms 
radii of the most and least localized Gaussians. The axes show distances in fm. 
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